arXiv: 1505.06680vl [physics.plasm-ph] 25 May 2015 


Optimized multibeam configuration for observation of QED cascades 


E.G. Gelfer, 1 A.A. Mironov, 1 A.M. Fedotov, 1 V.F. Bashmakov, 2,3 
E.N. Nerush, 2,3 I.Yu. Kostyukov, 2,3 and N.B. Narozhny 1 

1 National Research Nuclear University “MEPhI” (Moscow Engineering Physics Institute), 115409 Moscow, Russia 
2 Institute of Applied Physics, Russian Academy of Sciences, 603950 Nizhny Novgorod, Russia 
3 Lobachevsky National Research University of Nizhni Novgorod, 603950 Nizhny Novgorod, Russia 

QED cascades in intense electromagnetic field can occur if the dynamical quantum parameter \ 
of a seed electron, which in Compton units coincides with the electron proper acceleration, attains 
the order of unity. We derive general expression for x of an initially slow electron in an arbitrary 
electromagnetic field for a time range t -C 1/w, where u is the field carrier frequency. Using this 
formula, we consider a special field configuration of multiple colliding focused laser beams and 
optimize it to provide cascade development at laser power below 10 PW and intensity of the order 
of 10 23 W/cm 2 . Such parameters of the beams will be obtained with a new generation of laser 
facilities, particularly the ELI Beamlines, in the coming years. 
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I. INTRODUCTION 


During the last two decades one could observe a 
remarkable progress in laser technology. Since the in¬ 
vention of the CPA technique [l] in 1985, the max¬ 
imum laser intensity has been increased by seven 
orders of magnitude reaching nowadays the value 
2 • 10 22 W/cm 2 0- Moreover, the announced new gen¬ 
eration of 10 PW laser facilities like ELI Beamlines, 
which is now under construction in Prague, Czech Re¬ 
public 0 , or VULCAN in UK 0 , will provide inten¬ 
sity above 10 23 W/cm 2 . Therefore physics in the pres¬ 
ence of such ultra-strong fields is now coming upon a 
focus of attention. One of the most striking effects in 
this area is electron-positron pair creation from vac¬ 
uum M. However, according to recent calculations 
0, spontaneous pair creation from vacuum would re¬ 
quire an optical laser of the power above 500 PW, still 
far from the experimental possibilities of the nearest 
future. 

Another important example of ultra-strong field 
physics phenomena are quantum electrodynamical 
(QED) cascades [ToI 4T5| . In the present paper we con¬ 
sider cascades seeded by electrons, that are initially at 
rest at the focus of intense laser field 1 . They are ac¬ 
celerated by laser field emitting hard photons, which 
in turn create further electron-positron pairs. The 
created electrons and positrons are again accelerated 
by the field thus emitting hard photons, which again 
create new pairs and so on. Such an avalanche-like 
process finishes by escape of the particles from the 
focal region. 

Observation of QED cascade production in a labo¬ 
ratory is important for both fundamental physics and 
testing the QED principles. However, numerical sim¬ 
ulations performed in the previous works fioi . I© ] for 
some particular models of laser pulses have indicated 


that the threshold intensity required for cascade de¬ 
velopment would be of the order or above 10 24 W/cm 2 . 
In the latter case this may still exceed the capabilities 
of such facilities as e.g. the ELI Beamlines. Our goal 
is to optimize the experimental setup, in particular the 
field distribution at the focus, thus suggesting config¬ 
uration that decreases the threshold intensity required 
for QED cascade production. 

Consider the QED field strength parameter y [l8[ 
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where m and e are the electron mass and charge, e 
is the particle energy and Fj_ is the component of 
the Lorentz force per unit charge F = E + [v x B] 
transverse with respect to the particle velocity v. In 
the framework of a locally constant field approxima¬ 
tion 0,11m the parameter © evolves in between the 
hard photon emission and pair creation events accord¬ 
ing to the classical equations of motion and its instant 
value defines the probabilities of these processes fl8j 
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where W e and W 1 are the probabilities of pair creation 
and hard photon emission, respectively, e ej7 is the en¬ 
ergy of radiating electron or pair creating photon and 
a = e 2 //ic is the fine structure constant. Since y 7 is 
less than \ e of a radiating electron and pair creation 
is exponentially suppressed at y 7 < 1, we have to 
find such a field configuration that provides the most 
rapid growth of x 2 ■ Then we calculate numerically 
the threshold intensity for QED cascade development 
in such an optimal field configuration. 


1 Note that it is also possible to initiate QED cascades in a laser 
field by seed hard photons or high-energy electron beam [ifi j, 

the latter similar to the SLAC experiment D3. 2 Herein and onwards we abbreviate \e back to x • 
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II. GROWTH OF PARAMETER x 


Assume that a seeding particle is placed initially at 
rest at the focus of the laser field. In order to optimize 
field configuration for cascade production, we derive 
first a formula for the time dependence y(t) near the 
focus, i.e. in a time range t <C 1/w, where oj is the 
field carrier frequency. In this case the electric and 
magnetic fields can be expanded into Taylor series 

E(i) «E o + E ifc (0)a: fc , B(z) » B, fc (0)x* 

with respect to time and coordinates x k = {ct,x}. 
Since a standing wave composed of two identical coun- 
terpropagating laser beams produces cascades more 
efficiently than a single beam of the same total power 
0 , we assume in advance that the electric field at the 
center of the focal region reaches maximum and the 
magnetic field vanishes. Besides, in what follows we 
assume that the laser field is so strong that particles 
become ultrarelativistic soon after start of accelera¬ 
tion, eEt acc /mc 1, where t acc is the acceleration 
time, but still remains much weaker than the QED 
critical field, E <C E cr = m 2 c 3 /eh = 1.3 • 10 16 V/cm. 

Let us find the particle momentum satisfying the 
equation of motion p = eF in the form p = p (0 ) + 
pW + ..., where p( fe ) = 0{{ut) k ). For t m/eE 0 
3 the solution of equations of motion, accurate to the 
zeroth order in uit, reads 


p (0) = eE 0 t, e<°> = eEot, r<°> = v (0 >t = 

Eq 

Using expressions for the fields within the first order 
in wt, E(x) = Eo + E 't, B(x) = B 't, where 

E' = {|+ v(°) V J E(0), B' = j|- + v(°) v| B(0), 

( 2 ) 

we can solve the first order equation of motion 
dp^/dt = eFU) = ef^t, where 


fl 1 ) = E' + [ v (°) X B'], (3) 


and obtain for the particle velocity v = p/e, 
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Since the magnitude of the transverse component of 
the Lorenz force is equal to Fj_ = |[F x v]| = ef±\/ 2, 
where is the component of © transverse to the 
initial electric field Eq, we arrive at the formula 
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The result © is in agreement with special cases of 
rotating purely electric field and a standing plane wave 
obtained previously in Refs. El and 0. 

Hence, in order to optimize field configuration for 
cascade production we must increase the coefficient at 
t 2 in the formula © as much as possible. 


3 We use the natural units h = c = 1 onwards. 


III. COLLIDING BEAMS 



FIG. 1: (Color online). Allocation of n = 8 tightly focused 
beams colliding in a plane. 


A method for reducing the total threshold power 
for spontaneous pair creation from vacuum via multi¬ 
beam technology was suggested in Ref. where it 
was clearly demonstrated how the key parameter of 
that problem, the electric field strength in the focus, 
can be increased considerably by splitting a laser beam 
into several pairs of linearly polarized colliding beams 
of the same total power. However, in our case of cas¬ 
cade production linear polarization is expected to be 
unfavorable 00- Therefore we modify that scheme 
not restricting polarization of individual beams a pri¬ 
ory. 

Let us consider an even number n > 2 of focused 
beams propagating in xy plane at angles 


<Pj = 27r (j - l)Ab j = l--.n 

with the x axis (see Fig. © and optimize their po¬ 
larization for cascade production. For not too strong 
focused beams their field near the focus can be well 
approximated by plane waves 

A j =-iA 0 e j e- i ^ t - rAi l (6) 


Here e 3 = ej — ie J 2 , ej-Le^ are the individual polariza¬ 
tion vectors, which we do not require in advance to be 
normalized, and h 3 Ee 3 are the unit vectors directed 
along the focal axis of the beams. In case all the 
beams are focused weakly, the power of each of the 
beams can be estimated as Pj = Al\e 3 \ 2 />t 2 , where 
x = ^J&k/uj 2 Sx is a geometric factor depending on 
the beam shape, S± uj~ 2 is the focal cross section 
area of the beam. The same kind of formula remains 
true for a tight focusing case as well, but the geometric 
factor in such a case is modified (see below). Below, 
we prefer to fix the total power P = Pj of the 

beams, in terms of which 
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Due to transversality of the plane waves ©, all the 
terms containing space derivatives in © vanish within 
this approximation, so that B' = — *uB(0) = 0. Sub¬ 
stituting E = ^^dAj/dt into © and taking into 

j 







3 


account 0, for n colliding beams we obtain 
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where the factor 
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totally determines dependence of % on the number of 
beams and their polarization. 

Evidently, in order to maximize the expression @ 
all the vectors e{ should be parallel each other. Since 
the cry-plane had been chosen as the plane of beam 
propagation, we assume thus that e\ are all directed 
along the 2 -axis, = gji , where z is the correspond¬ 
ing versor. Also, since the electric field in the beam is 
transverse to the wave propagation direction rP 4 , we 
can assume that e 2 = —hj sin^-x + hj cos^y. 

Due to the symmetry of the problem the polariza¬ 
tion vector of the total electric field can lie in any 
plane parallel to 2 -axis. Let it be the yz- plane. Then 
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where, to clarify further manipulations, in the last rep¬ 
resentation we have introduced four formal n-vectors 
G = {gj}, H = {hj}, a = (cos ipj} and (3, the latter 
with all n components equal to 1. 

Now it follows from the chain of inequalities 5 
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that T n reaches its maximum value when G || (3, H || 
a and ||G|| = ||H||. For n > 2 these conditions are 
equivalent to 

g 3 = gi, £i = ngiz, hj = hi cos ipj, (12) 

n 

£2 = hiy ^2 cos 2 <fij = —/iiy. (13) 

3=1 

In particular, 

ei : e 2 = V2 : 1, (14) 

i.e. the optimal polarization of the total field is ellipti¬ 
cal with the semi-major to semi-minor axis ratio equal 


4 Note that since the magnetic field in the focus is assumed to 

be equal to zero, we must impose a constraint , 

or just hj = —h n / 2 +j, for all j = 1... n/ 2. 

5 Here we apply by turn the inequality of arithmetic and geo¬ 
metric means and then the Cauchy-Bunyakovsky inequality. 


to y/2 : 1. The polarization vectors of the individual 
beams in such an optimal case are 

€j = {iV% sin (fj cos ipj , —iy/2 cos 2 (pj , 1}. (15) 

In case n = 2 the factor 2 does not appear in denom¬ 
inator of Eq. da, hence the optimal polarization is 
circular, e\ : e% = 1 : 1. 


IV. NUMERICAL RESULTS 


TABLE I: Power and intensity thresholds for generation 
of QED cascades in configurations with n colliding beams 
(w = leV, £ = 8). 


n 

2 

4 

6 

8 

Pth, PW 

23 

14 

10 

7.9 

Ith, 10 23 W/cm 2 

17 

10 

7.3 

5.6 



FIG. 2: (Color online) Number of electrons in a cascade 
versus time for optimal (solid) and circular (dashed) po¬ 
larization of the total field (n = 8, £ = 8, uj — leV, 
P = 7.9 • 10 15 W, T is the laser period). 


The threshold power for cascade generation can be 
found by numerical simulations. In fact, QED cas¬ 
cades do not possess a natural sharp threshold, so 
that we define it formally as creation of just a single 
pair on average by each seed electron during passage 
of laser pulses through the focus. As a more realis¬ 
tic model we use a model of the focused laser beam 
suggested in Refs. @, l20| 
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FIG. 3: (Color online), (a) number of electrons at a distance smaller than A from the origin; (b), (c) - dimensionless 
linear densities A of created positrons along y- and z-axes. The parameters are the same as in Fig. [2] 



FIG. 4: (Color online). The mean value of \ for electrons 
located at a distance smaller than A from the origin. The 
parameters are the same as in Fig. [2] 

which is an approximate solution of Maxwell equa¬ 
tions near the focal region if £ 1. Here, and 

are the components of the radius vector, parallel and 
perpendicular to n-h The parameter £ = uL/ 2, where 
L is the diffraction length, controls the degree of beam 
focusing; rV, e•?' and e are all defined after Eq. (j6]); P 
is the total power of all colliding beams. 

It follows from (JSJ that for given total power, the ge¬ 
ometric factor x (i.e., focusing) should be maximized. 
Note that in modern laser facilities the beam can be 
focused so tightly, that a half of the beam power is 
passing through a circle of the wavelength diameter in 
the focal plane @ . Computing numerically the power 
of a beam over a circle of the diameter A = 2i t/uj 
around the focus, we obtain, that the condition 

r \/2 

T Re /*/ pdp ft > ■ [E, X B *]| rjj _ 0 = h 

0 0 

where Pj is the power of the j-th beam, is satisfied if 
the focusing parameter £ ~ 8. Here p and (j) are the 
polar coordinates in the plane = 0. 

It is also worth taking into account that there exist 
bounds for a number of beams that can be allocated 
concurrently in a plane without overlapping @,ET|. At 
far distance from the focal region the field of each of 
the beams is almost confined to a cone with aperture 
0, which is simply related to our focusing parameter 
£. Indeed, for |r| £/w the energy of the field is 

proportional to e 4r - L/r ii. Therefore we can estimate 


the aperture by 9 ~ 2 arctan 1/-y/£. For £ = 8, 9 = 
O.2107T > 7t/ 5, hence the number of beams that can 
be allocated without overlapping in the :cy-plane must 
be less than 10. Hence, for £ = 8 the maximal even 
number of beams is 8 as in Fig. [I] 

For numerical simulations we have used two inde¬ 
pendent codes. The first one [l6| is a realization 
of Monte-Carlo (MC) algorythm, in which we as¬ 
sume that electrons and positrons are moving clas¬ 
sically between the emission events, which are imple¬ 
mented by the corresponding event generators similar 
to Refs. [12, Hj]. The second one is a more general 
3D Particle-in-Cell-Monte-Carlo (PIC-MC) code for 
simulation of laser-plasma interactions with quantum 
radiation effects [l3j, [01 (221 • Actually, in our problem 
the electron-positron plasma remains dilute and the 
results obtained with both codes are in good agree¬ 
ment. In this paper we present the results obtained 
with the first (pure MC) code. 

Since our field model m corresponds to a field 
unbounded in time, we take into account finiteness of 
pulse duration artificially by limiting the simulation 
time. The simulation time was chosen to be 5 laser 
periods T as is typical for modern laser facilities Q. 
In simulations, the seed particles were represented by 
a bunch of Nq = 10 3 electrons placed initially at rest 
near the center r = 0 of the focus. The laser car¬ 
rier wavelength was taken to be A = 1.24/xm, which 
corresponded to laser photon energy leV and the fo¬ 
cusing parameter was chosen as £ = 8 for the reasons 
discussed above. 

The resulting power and intensity thresholds ob¬ 
tained by numerical simulation of QED cascade de¬ 
velopment in the cases of n = 2, 4, 6, and 8 beams 
with optimal polarization (1151) . colliding in a plane are 
presented in Table [U One can observe that both the 
laser power and intensity required for cascade produc¬ 
tion decrease significantly due to a multibeam setup. 
Cascade dynamics for the cases of optimal elliptic po¬ 
larization and circular polarization of the total field 
for the case of n = 8 beams is displayed on Fig. [2] 
The figure clearly demonstrates the advantage of the 
suggested optimal polarization. 

Besides, it is clear from Fig. [2] that in the case of 
optimal polarization cascade development is splitted 
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into two stages, which differ by the slope of the curve. 
At the first stage the seeding particles are accelerated 
by the field, create pairs and escape the focal region. 
Almost all created particles also leave the focal region 
with a speed close to the speed of light, see Fig. [3] 
Still, a small fraction of slow particles remains in the 
focal region being accelerated by the field. At time 
t ~ 2.25 T the mean value of the parameter \ for par¬ 
ticles at the focus becomes large enough to initiate the 
second stage of cascade development, see Fig. 2) 

V. CONCLUSION 

We have demonstrated that QED cascades can arise 
in laser-target interaction for laser power below 10 PW 
and intensity below 10 24 W/cm 2 . This is realizable 
with a multibeam setup, i.e. splitting of the origi¬ 
nal laser pulses into several pairs of colliding beam. 
In order to determine how to arrange the beams, we 
have derived an analytical formula for growth of the 
parameter x(t), which determines possibility of QED 
cascade generation. This allowed us to identify the op¬ 
timal beam polarization with the most rapid growth 
of x(t) at the focus. For such an optimal configu¬ 
ration of n beams and fixed total laser power, x is 
proportional to n. However, for our extreme choice 


of the focusing parameter, which corresponds to the 
state-of-the-art experimental capabilities, the maxi¬ 
mal number of beams that can be brought together 
in a planar configuration of Fig. [T] is n = 8. Thus 
we have next computed the QED cascade thresholds 
for the particular case of n = 8 optimally polarized 
beams. According to our simulations, in this case the 
threshold total power P t h = 7.9 PW and the thresh¬ 
old intensity I t h = 5.6 • 10 23 W/cm 2 . Such parameters 
had been previously announced as being achievable 
with the new generation of laser facilities SB in the 
nearest future. 
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